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In this paper we present a variant of the Calculus of Looping Sequences (CLS for short) with global 
and local rewrite rules. While global rules, as in CLS, are applied anywhere in a given term, local 
rules can only be applied in the compartment on which they are defined. Local rules are dynamic: 
they can be added, moved and erased. We enrich the new calculus with a parallel semantics where 
a reduction step is lead by any number of global and local rules that could be performed in parallel. 
A type system is developed to enforce the property that a compartment must contain only local rules 
with specific features. As a running example we model some interactions happening in a cell starting 
from its nucleus and moving towards its mitochondria. 

1 Introduction 

The Calculus of Looping Sequences (CLS for short) @|4l[T5][6), 1S a formalism for describing biological 
systems and their evolution. CLS is based on term rewriting with a set of predefined rules modelling the 
activities one would like to describe. CLS terms are constructed by starting from a set of basic con- 
stituent elements which are composed with operators of sequencing, looping, containment and parallel 
composition. Sequences may represent DNA fragments and proteins, looping sequences may represent 
membranes, parallel composition may represent juxtaposition of elements and populations of chemical 
species. 

The model has been extended with several features such as bisimulations El 13, combining the sim- 
plicity of notation of rewrite systems with the advantage of a form of compositionality. In H] 13 a type 
system was defined to ensure the well-formedness of links between protein sites within the Linked Cal- 
culus of Looping Sequences (see 0). In lfT3l we denned a type system to guarantee the soundness of 
reduction rules with respect to the requirement of certain elements, and the repellency of others. 

In this paper we present a variant of CLS with global and local rewrite rules (CLSLR, for short). 
Global rules are applied anywhere in a given term wherever their patterns match the portion of the 
system under investigation, local rules can only be applied in the compartment in which they are denned. 
Terms written in CLSLR are thus syntactically extended to contain explicit local rules within the term, on 
different compartments. Local rules can be created, moved between different compartments and deleted. 
We feel that having a calculus in which we can model the dynamic evolution of the rules describing the 
system results in a more natural and direct way to study emerging properties of complex systems. As it 
happens in nature, where data and programs are encoded in the same kind of molecular structures, we 
insert rewrite rules within the terms modelling the system under investigation. 

In CLSLR we also enrich CLS with a parallel semantics in which we define a reduction step lead by 
any number of global and local rules that could be performed in parallel. 
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Figure 1 : From the cell nucleus to mitochondria. 

Since in this framework the focus is put on local rules, we define a set of features that can be asso- 
ciated to each local rule. Features may define general properties of rewrite rules or properties which are 
strictly related to the model under investigation. We define a membrane type for the compartments of 
our model and develop a type systems enforcing the property that a compartment must contain only local 
rules with specific features. 

Thus, the main features of CLSLR are: 

• different compartments can evolve according to different local rules; 

• the set of global rules is fixed; 

• local rules are dynamic: they can be added, moved and erased according to both global and local 
rules; 

• a parallel reduction step permits the application of several global and local rules; 

• compartments are enforced to contain only rules with specific features. 

As a running case study, emphasising the peculiarities of the calculus, we consider some mitochon- 
drial activity underlining the form of symbiosis between a cell and its mitochondria (see lfl4lD . Mitochon- 
dria are membrane-enclosed organelle found in eukaryotic cells that generate most of the cell's energy 
supply in the form of adenosine triphosphate (ATP). A mitochondrion is formed by two membranes, the 
outer and the inner membrane, having different properties and proteins on their surfaces. Both mem- 
branes have receptors to mediate the entrance of molecules. In Figure Q] we show the expression of a 
gene (encoded in the DNA within the nucleus of the cell) destined to be translated into a protein that will 
be catch by mitochondria and will then catalyse the production of ATP. In particular, we will model the 
following steps: (1) genes within the nucleus' DNA are transcribed into mRNA, (2) mRNA moves from 
the nucleus to the cell's cytoplasm, (3) where it is translated into the protein. 

The vast majority of proteins destined for the mitochondria are encoded in the nucleus of the cell and 
synthesised in the cytoplasm. These are tagged by a signal sequence, which is recognised by a receptor 
protein in the Transporter Outer Membrane complex (TOM). The signal sequence and adjacent portions 
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of the polypeptide chain are inserted in the intermembranous space through the TOM complex, then in 
the mitochondrion internal space through a Transporter Inner Membrane complex (TIM). According to 
this description, we model the following final steps: (4) protein is detected by TOM and brought within 
the intermembranous space, (5) then, through TIM, in the mitochondrion's inside, (6) where it catalyses 
the production of ATP, (7) that exits the inner, (8) and the outer mitochondrial membranes towards the 
cell's cytoplasm. 

2 The calculus 

In this Section we present the Calculus of Looping Sequences with Local Rules (CLSLR). 
2.1 Syntax of CLSLR 

We assume a possibly infinite alphabet S of symbols ranged over by a,b,c, . . ., a set of element variables 
3£ ranged over by x,y,z, . . ., a set of sequence variables ranged over by x,y,z, . . ., and a set of term 
variables 3"Y ranged over by X,Y,Z, . . .. All these sets are possibly infinite and pairwise disjoint. We 
denote by f the set of all variables, Y = JT U ,9"Y U STY , and with % a generic variable of Y. Hence, 
a pattern is a term that may include variables. 

Definition 2.1. [Patterns] Patterns P, sequence patterns SP and local rules R o/CLS are given by the 
following grammar: 



p 


::= (SP)°\P P\P 


L 


SP 


::= e a x SP 


SP | x 


L 


::= X SP L L 


R 


R 


::= L^L \ L^ SP ^L^ SP 


| l Isp^ l isp 



where a is a generic element of S, and X,x and x are generic elements of £7"f and S£ ', respec- 
tively. Sequence patterns SP defines patterns for sequences of elements, (SP) denotes a closed (looping) 
sequence which may contain other patterns through the J operator, | is used to denote the parallel com- 
position (juxtaposition) of patterns, and R denotes the syntax of local rules that may either exit (L^ SP ) or 
enter (L^ SP ) a closed sequence SP. We denote with & the infinite set of patterns. 

A local rule L\ ^tLi is well formed ifL\ ^ £ and Var(L,2) C Var(L\), where Var(P) denotes set of vari- 
ables appearing in P. 

A local rule l\ sp ' ^l\ SPi or l\ SPi ^L^ Pl is well formed ifU # £, Var(L 2 ) C Var(L x ) and Var(SP 2 ) C 
Var(SP { ). 

Terms are patterns containing variables only inside local rules. Sequences are closed sequence pat- 
terns. We denote with S' the infinite set of terms, ranged over by T, and with S" the infinite set of 
sequences, ranged over by S. 

An instantiation is a partial function a : (3"f ->• 3") U (S^Y ->• S?) U ( X ->■ S). Given P£#, 
with Pa we denote the term obtained by replacing each occurrence of each variable appearing in 

P with the corresponding term o(%), but for local rules, which are left unchanged by instantiations, i.e., 
Ra = R for all 7? and a. With £ we denote the set of all the possible instantiations. 

Definition 2.2 (Structural Congruence). The structural congruence relations =$ and =# and =p are the 

least congruence relations on sequence patterns, local rules and on patterns, respectively, satisfying the 
rules shown in Figure |2] 
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SI\ ■ (SP 2 • SP 3 ) = s (SPi -SP 2 )- SP 3 SP e= s e SP= s SP 
SPi = s SP 2 implies SPi = P SP 2 and (SP])° \P= P (SP 2 f}P 

L\ =p Lj and L 2 =pL' 2 imply L\ ^L 2 =p L\ h^-L' 2 
Li= P L\ and L 2 = P L' 2 and SPi= P SP'i and SP 2 = P SP' 2 
imply Lj SP ' ^Lf 2 =* Z/f 1 H+Z/f" 2 and l\ sp > ^L^ = r L'f P '< ^L'f P ' 2 

7?i =rR 2 implies Pi =p R 2 
P l \P 2 = P P 2 \P P 1 \(P 2 \P 3 )= P (P 1 \P 2 )\P 3 P\e= P P 
(e)° J £ = P £ (5Pi • SP 2 )° J P = P (SP 2 • SPi) j P 

Figure 2: Structural Congruence 

Structural congruence rules the state the associativity of • and | , the commutativity of the latter and 
the neutral role of £. Moreover, axiom (SP\ ■ SP 2 ) J P =p (SP 2 ■ SP\ )° J P says that looping sequences 
can rotate. In the following, for simplicity, we will use = in place of =p. 

2.2 (Parallel) Operational Semantics 

In order to define a reduction step in which (possibly) more than one rule is applied, following iflOl . we 
first define the application of a single rule, either global or local. We resort to the standard notion of 
evaluation contexts. 

Definition 2.3 (Contexts). Evaluation Contexts E are defined as: 

E::=D \ E\T j T\E \ (S)°\E 

where T € ST and S G 5^ . The context □ is called the empty context. We denote with E the infinite set 
of evaluation contexts. 

By definition, every evaluation context contains a single hole □. Let us assume E G E, with E[T] we 
denote the term obtained by replacing □ with T in E. The structural equivalence is extended to contexts 
in the natural way (by considering □ as a new and unique symbol of the alphabet $"). Note that the shape 
of evaluation contexts does not permit to have holes in sequences. A rewrite rule introducing a parallel 
composition on the right hand side (as a i->- b \ c) applied to an element of a sequence (e.g., m-a-m) would 
result into a syntactically incorrect term (in this case m ■ (b \ c) ■ m). 

To enforce the fact that local and global rule applications can be done in parallel, we underline those 
subterms that are produced by the application of the rule involved. Terms matching the left4iand-side 
of a (local or global) rule must not have any underlined subterm. Underlined terms are only used for 
bookkeeping in the definition of rule application. 

Let 2?_ denote the set of terms in which some subterms can be underlined. The erasing mapping 
T] : 3F_ i-> Sf erases all underlining obtaining a term generated by the grammar of Definition 12. II 

We first define the application of global or local rules to terms in 3_ which produces terms in £F . 

Global rules are of the shape Pi i-> P 2 . They can be applied to terms only if they occur in a legal 
evaluation context. 

Local rules are inside compartments, and can be applied only if a term matching the left4iand-side 
of the rule occurs in the same compartment. 

Notice that global rules have patterns P on both the left and the right4iand side of the rules, whereas 
local rules have the less general L patterns. In particular, L patterns do not contain compartments, and 
therefore cannot change the nesting structure of the compartments of a term. 
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Definition 2.4 (Rule Application). Given a finite set of global rules Si, the rule application —> is the least 
relation closed with respect to = defined by: 

Pi^P 2 e& ctg! P\0^e 

(GRT) 



E[Pio] ^E[P 2 a] 
a G £ L\G ^ e 



(LR) 



E[L\ ^L 2 \ L x o | T] -+E[L\ ^ L 2 \ L 2 o_ \ T] 
a G I LiO ^ £ 

£[(5ia) 1 (Zf ^Lf 2 | L x a \ T)]^E\L 2 a \ (S 2 o)° \ {If 1 ^ Lf 2 \ T)} (LR "° UT) 

E\L\ Sl ^Lf 2 | Uo | (5ig)° | T]^E\L\ Sl ^Lf 2 | (5 2 g)° | (T \ L 2 o)\ 

With rule (lr-Out) a term or a rule Lid exits from the membrane in which it is contained only if the 
membrane has as loop the sequence Si a: when outside, L\G is transformed into L 2 a, and the sequence 
Si a into S 2 o. (lr-in) is similar, but it moves terms or rules into local membranes. Local rules do not 
permit to move, create or delete membranes: only global rules can do that. 

Observe that local rules can be dynamically added and deleted both by global and local rules. A 
global rule which adds the local rule R is of the shape P h-> R and a global rule which erases the same 
local rule is of the shape R h-> P. A local rule which adds the local rule R is of the shape Li-> R and a 
local rule which erases the same local rule is of the shape R i— > L. Moreover local rules can add local 
rules in compartments separated by just one membrane, since they can be of the shapes lJ s i-> R^ s or 

A reduction step of the parallel semantics =>, starting from a term in 3F applies any number of 
global or local rules that could be performed in parallel, producing a final term in & (with no underlined 
subterms). 

Definition 2.5 (Parallel Reduction). The reduction between terms in 2? is defined by: 

T = r () -> Ti ->•■•-> r w+ i n>0 T' = rj{T n+1 ) 
T =>T' 

To justify the definition of the reduction we have to show that the order in which the local or 
global rules are applied is not important. The notion of multi-hole context, i.e., of term where some 
disjoint subterms are replaced by holes is handy. More precisely the syntax of multi-hole contexts is: 

C::=D | T \ C\C \ (□)°JC | (S)° \C 

We can show that in a parallel reduction only disjoint subterms can change. Rules (GRT) and (LR) 
modify just one subterm. Rules (lr-out) and (lr-in) modify three subterms, i.e., a membrane and a term 
exiting or entering the membrane, using e for the missing term. 

Theorem 2.6. If T ==?■ T', then there is a multi-hole context C[ ] . . . [ ] and terms T\ , . . . , T n , . . . , T' n 
such that T = C[Ti] . . . [T n ], T 1 = C[T[] . . . [T^\ and for all 1 < i < n: 

. either C[T{] . . . [7}] . . . [T*\ -> C[T{] . . . [T{] . . . [T*] 
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• or there are i\ , i 2 , h such that i € {i\ , i 2 , i$ } and 

C[T*] . . . [T h ] | ([T h ]f\ [T h ] . . . [T n *] -> C[T*} . . . [7/J | ([7£])° j [TQ . . . [T*] 

where TJ can be either Tj (subterm ofT)or Tj (subterm ofT'). 

Proof. If T T', then for some Uq,--., U m+ \ we get T = Uq — > • • • — > U m+ \ where m > and T' = 
T^]{U m +i)- We show by induction on h < m and by cases on the last applied reduction rule that 

. either U h = C[T*\ . . . [7}] . . . [T*] and U h+1 = C[T*\ . . . [T[] . . . [T*\ 

• or there are i\ , ii , z'3 such that i £ { ?i , i 2 , 13 } and 

U h = C[T*] . . . [T h ] I ([T i2 ]f j [7} 3 ] . . . [T n *] and U h+l = C[T*] . . . [7/J | ([TQ) J [TV] . . . [T*] 

where " * " can be either " " or " ' " and all terms with ' are either underlined or e. 
If the last applied rule is 

a G £ L\G ^ e 
E[Li ^ L 2 I LiO I V] -> E[Li h-> L 2 \ UP_ \ V] 

then U h = E[L l ^L 1 \Lio\ V] and U h+1 = E[L X ^L 2 \ L2O \V]. By induction U h = C[T*\ . . . [T*]. 
Since L\0 is a subterm of T and L 2 a is a subterm of f/ m +i there must be an index i such that Ti = L\G 
and T{ = L 2 a . 

If the last applied rule is 

ael La^e L x c e ST S\G g 2? 
E\L\ Sl ^ Lf 2 I L t a I (fra) | V] -> £[Zf 1 ^ | fea)° | (V | 

then % = E[L\ S> ^ Lf 2 | Li<7 | (Sid) j V] and C^+i = E[l\ S 1 ^ Lf 2 | (S 2 g)° | (V | Z^cr )]. By in- 
duction % = C[T*] . . . [T*]. Notice that LiO,SiO are subterms of T and L 2 a , S 2 a are subterm of U m+ i. 
Moreover L\0, S\0 and e in L are replaced by e, S 2 G and L 2 o in U m+ i, respectively Therefore there 
must be indexes 11,12, 13 such that 7}, = L\0, T[ x = e, T 2 = S\G, T- 2 = S 2 a , 7] 3 = e, T^ = L 2 a . 

□ 

Example 2.7. [Mitochondria Running Example: Syntax and Reductions] A CLSLR term representing 
the mitochondria evolution inside the cell's activity discussed in the introduction could be: 

CELL = (cell) J ( NUCLEUS \ MITOCH \ ... \ MITOCH \ 
mRNA I—?- protein \ 
protein^- 70 ™ H> protein^- 7 '"" ) 

A cell is composed by its membrane (here just represented by the element cell) and its content ( in this 
case, the nucleus, a certain number of mitochondria and a few rules modelling the activity of interest). 
In particular, the two rules above model the steps (3) and (4), respectively, of the example schematised 
in the introduction. 

Assuming that DNA is the sequence of genes representing the cell's DNA, and g is the particular 
gene (contained in DNA) codifying the protein, we define the nucleus of the cell with the CLSLR term: 

NUCLEUS = (nucleus) j ( D N A | 

x ■ g ■ (x-g-y\ mRNA) \ 
mRNA^ nucleus \— > mRNA^ nucleus ) 
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(cell) ] ((nucleus) 



(cellf_ 
(cell) 



(mRNA | 
(protein 



(cell) ] ((Torn) 
(cell) \ ((Tom) 
(Tom) 
(cell) ] ((Tom) 

(Tom) 
(cell) J ((Tom) 
(Tom) 
(cell) ] (ATP\ .. 

(Tom) 



°] (mRNA | ... | mRNA | ...) | ...) 
.. | mRNA | ...) 
. . . | protein | .. .) 

(protein \ ...) \ ... \ (Tom) J (protein \ ...) 
((Tim) J (Mit A (Mif A | A77>) | . . .) | . . .) | 
((rim) J (Mit A ^(Mit A \ATP) | ...) | ...) | . 

((rfei)°j(A7P|...)|...)l 

((r/mf j(AJ J P| ...)[...)[...) 
(Arr I (Tim) J (A7P | ...) | ...) | 
(ATP | (rim) J (ATP \ . . .) j . . .) | . . .) 

| Arr | (rom)°j (Arr | (rim) J (Arr | . . .) 

(AIP| (rim) J (ATP | ...) | ...) | ...) 



...) 



Figure 3: Mitochondria evolution 



Note that the first of the two rules above models step (1) of our example (DNA transcription into mRNA ), 
the second one (mRNA exits the nucleus) models step (2). 

The mitochondria of our model are composed of a membrane, on which we point out the Tom com- 
plex, containing an inner membrane ( INN_MITOCH) and a couple of rules: 

MITOCH= (Tom) ]( INN_MITOCH\ 

protein^ Tim H- (Mit A i-> (Mit A \ ATP)f Tim \ 
ATP^^ATP^ ) 

where we denote with the element Mit A a mitochondrial factor inside the inner membrane ( activated by 
our protein ), necessary to produce ATP. In particular, the protein, when in the intermembranous space, is 
moved through Tim inside the inner mitochondrial space (step (5) of our example) and then transformed 
into the newly generated rule Mit A i— >■ (Mit A \ ATP) which will lead the production of ATP (step (6) of our 
example). 

Finally, in INN_MITOCHwe point out the Tim complex: 

INN_MITOCH= (Tim) J ( Mit A \ 

ATP^^ATP^ ) 

Both in MITOCH and INNJ4IT0CH we have the rules needed to transport the ATP towards the cell's 
cytoplasm (steps (7) and (8) of the example). 

A possible (parallel) reduction of this term, when g initially produces a certain number of mRNA is 
(by focusing only on the more interesting changes) shown in Figure\3\ The ATP produced in the last but 
two reductions in INN_MITOCH moves to MITOCH in the last but one reduction and new ATP is produced 
in INN_MITOCK In the last reduction, the firstly generated ATP moves to the cell, the secondly generated 
ATP moves to MITOCH and new ATP is produced in INNJtlTOCH 



3 Types 

In this section we inttoduce a type system that enforces the fact that compartments must contain rules 
having specific features. E.g., in [161 the following rule features for L\ H> L2 are suggested: 
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• the rule is deleting if Vars{L\) D Vars(L 2 ) (denoted by d); 

• the rule is replicating if some variable in L2 occurs twice (denoted by r); 

• the rule is splitting if L\ has a subterm containing two different variables (denoted by s); 

• the rule is equating if some variable in L\ occurs twice (denoted by e). 

This kind of features reflects a structure of rewrite features which could be common for rewrite systems 
in general. Other, model-dependent, features could be defined to reflect peculiarities and properties 
of the particular model under investigation. The features of the rules allowed in a compartment are 
controlled by the wrapping sequence of the compartment. Our typing system and the consequent typed 
parallel reduction ensure that, in spite of the facts that reducing a term may move rules in and out of 
compartments, compartments always contain rules permitted by their wrapping sequence. In addition to 
the previous features of rules we say that: 

• the feature of rule L\ Sl 1— >• L^ 2 is that it is an out rule (denoted by o); 

• the feature of rule l\ S[ 1— >■ L^f 2 is that it is an in rule (denoted by i). 

To express the control of the wrapping sequence over the content of the compartment, we associate a 
subset cp of {d, r,s,e,o, i} to every element in $ . This is called a membrane type. We use A to denote a 
classification of elements. The type assignment in Figure 01 where a basis A assigns membrane types to 
element and sequence variables, defines the type of a sequence as the union of the membrane types of its 
elements. 





a : (p G A 


A h s £ : (TSeps) A,X (p\~ s X '• <P (TSvar) 


(TSelm) 


A h .5 a : (p 


A h 4 SP : «p A h s SP' : (p' 




(TSSEQ) 




Ah s SP SP' : <pU«p' 





Figure 4: Typing Rules for Membranes 

To define the type of patterns, that may contain parallel (composition) of rules, we consider: 

1. the features of the rules, contained in the pattern, and 

2. in case there are output rules the type of the rules that are emitted by these output rules. 

Therefore a pattern type, denoted by T, is a sequence of membrane types, i.e., T S {<?>}*■ With we 
denote the empty sequence. If the parallel composition of local rules R\ \ ■■■ \ R n , n > 0, has type (p ■ T, 
then (p is the union of the features of the rules Ri (1 < i < n), and z is the type of the parallel composition 
of rules in L' for those L' such that /?, = L^ SP \— > L'^ SP for some i, 1 < i < n (the type of the parallel 
composition of the rules that are emitted). If no rule is emitted, then t = 0. 
Union of pattern types, U, is defined inductively by: 

• 0Ut = tU0 = t, and 

• <pi • Ti U (p2 • T 2 = Ol U (p2) • (Tl U T 2 ). 

and containment, C, is defined by: 

• 0CT 
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• 9 • T C <p' • t' if <p C <p' and T C t'. 
The judgment A hp P : z, denned in Figure [51 asserts that the pattern P is well formed and has pattern 
type Z, assuming the basis A, which assigns membrane types to element and sequence variables and 
pattern types to term variables. The judgement A \- gr P\\-^P 2 : ok in last rule defines well-formedness of 
global rules. It is easy to verify that the typing rules in Figures |4] and [5J enjoy weakening, i.e., if A C A' 



A,X : z\- p X : z (Tvar) A \- p SP : (Tseq) 

A hp L 2 : T 

(TRLOC) 

A hp L\ H> L 2 : features(Li i— > L 2 ) U t 
A hp L 2 : t A h, Si : <pi A h, S 2 : (jh (piQ<(>2 



A hp Lf 1 i-> L| 52 : {o} • t 
A hp L 2 : <p • z A \- s 5i : <pi A h s . S 2 : <p 2 <p U <pi C ^ 



(TRlocOut) 



AhpLf 1 ^Lf 2 : {i}UT 



(TRlocIn) 



A hp Z 5 : t A hp P' : t' 
AhpP|P':TUT' 



(Tpar) 



A h 4 SP : (p A hp P : (p' ■ z' (p' C <p 
Ah p (SP)°J/ , :T / 

A hp A : Ti A hp P 2 : T 2 T 2 □ n 



(TCOMP) 



A h gr Pi^-P 2 : ok 



(TRglob) 



Figure 5: Typing Rules for Patterns and Global Rules 

then A h s SP : <p implies A' h 4 SP : <p, A h p P : t implies A' h p P : z, and A \- gr Pi i-» P 2 : ok implies 
A' h gr Pi^P 2 : ok. 

Rule (Tvar) asserts that a term variable is well typed when its pattern type is found in the basis. Rule 
(Tseq) asserts that, since a sequence does not contain rules, its pattern type is empty. Rule (TRloc) asserts 
that the type of a local rule R = L\ i-> L 2 is the union of the set of features of the rule R, denoted by 
features(P), and the pattern type of its right-hand-side L 2 . This is because once the rule is applied an 
instance of the pattern L 2 will substitute the instance of its left-hand-side L\. Rule (TRlocOut) checks 
that the features of rules permitted by the membrane S 2 include the one permitted by Si , so that if the 
compartment was well formed before applying L[ Sl (->■ L^ 2 , it will be well formed afterwards (when S 2 
replaces Si). Moreover, the pattern type of the rule is {o}, concatenated with the pattern type of L 2 , since 
L 2 is the pattern sent outside the compartment. Rule (TRlocIn) checks rule l\ S[ h-> df 1 . Since the pattern 
L 2 will get into a compartment with membrane Si, the membrane S 2 , that replaces Si, must permit all the 
features of rules that were permitted before, and moreover, it permits the features of the rules in L 2 . The 
type is {i} union the type of the patterns that are emitted by the out rules contained in L 2 . Rule (Tpar) 
enforces the fact that the patterns in parallel are both well formed and the final pattern type is the union 
of the two pattern types. Rule (Tcomp) checks that a compartment contain only rules whose features are 
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permitted by its wrapping sequence. The pattern type of the compartment is the type of the rules that are 
emitted. Finally, rule (trglob) says that the global rule Pi h-> P 2 is well formed in case the pattern P2 that 
will replace Pi has less features, so that it is permitted by all the compartments in which Pi is permitted. 

As we can see from rule (trloc) the type system is independent from the specific set of features 
considered. Any syntactic characterisation of rules could be considered for a feature. 

Based on the previous typing system we define a typed semantics, that preserves well-formedness of 
terms. Let an instantiation a agree with a basis A (notation a G Ea) if x : <p G A implies h s a(x) : cp, 
x : q> G A implies \- s a (5c) : (p, and X : z G A implies \- p o(X) : z. This is sound since the judgments 
\- s use assumptions on element and sequence variables, while the the judgments \- p use assumptions on 
term variables. 

Definition 3.1 (Typed Rule Application). Given a finite set of global rules M, the typed rule application 

— > T is the least relation closed with respect to = and satisfying the following rules: 

^ A = {Pj P 2 G M I A h gr Pi ^ P 2 : ok} 
Pi H- P 2 G a G I A Pl<7 # £ 

(T-GRT) 

£[Pia] ^ T E\P 2 a} 
(J€l A Lia^e 

(T-LR) 



£[Li i-^L 2 I L\a] -^jE[Li h>L 2 | L 2 a ] 



(T-LR-OUT) 



£[(5ig) d I (L] s 1 ^L^ 2 I L t g I P)1 ^ T P[L 2 g | (S 2 g)° | (Lf 1 ^ Lf 2 \ T)\ 

ctgIa L\G^e 

E\L\ S> ^pf 2 I Lio I (5ig)° I rj ^ t E\l\ s> ^L\ Sl I (S 2 g)° | (T \ L 2 o)] (T " LR " lN) 
Definition 3.2 (Typed Parallel Reduction). P/je reduction ==> T between term in is defined by: 

T = T ^jT l ^ J ---^ J T n+l n>0 T' = T](T n+1 ) 

T' 

The property enforced by the type system is that well-typed terms reduce to well-typed terms: the 
proof is the content of the Appendix. 

Theorem 3.3 (Subject Reduction). If\- p T:z and T =^y T', then h p T' : z' for some z 1 C Z. 

Example 3.4 (Mitochondria Running Example: Typing). Let the rules used in Example \2.7\ be labelled 
by: 

• R g = x ■ g ■ jTi-4 (5c-g-y\ mRNA), 

• R m = mRNA 1— > protein, 

• R ^ = protein* Tom 1— > protein^ 70 " 1 , 

\Tim 



• P mt = mRNA^ nudeus ^mRNA^ nudeus , 



• P,; = protein^ Tim ^Rj 



P fl = Mit A ^ (Mit A I ATP) 
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• =ATP^ x ^ATP tx . 

Let (p g = features(7?g), <p m = features(/? m ), (p a = features(/? u ). The term representing our model can be 
typed if A contains appropriate membrane types for the elements which occur in the membranes, i.e.: 

{cell : (p ce! i, nucleus : (p mc ieus,Tom : (p Tom ,Tim : (p T im} Q A 

where {i} U (p,„ C (p ce u, {o} U (p g C % uc i eus , {o, i} C (pjom, and {o} U (p a € <Pn'm- In this case the given 
parallel reduction is also a typed parallel reduction for this term. 
We can type the MITOCHwith the following derivations: 



A \~ p Mit A | ATP : 

(TRLOC) 

A \- p R a : (p a A h s Tim : (p Tim (p a C (p Tim 

(TRLOClN) 

AH p %:{i} Ah p R 4 :{o} 

■ ; ■ (TPAR) 

Ah p % \R cn : {i,o} 



A h p INNJUTOCH : A h p % | R a ^ : {i, o} 

(TPAR) 

A \- s Tom : <p rom A h p INNJIITOCH \ R^ \ R^ : {i , o} { i , o} C <p rom 

7Z 1 ! (TCOMP) 

A hp (Tom)° J (INN_MIT0CH\ % | /? at ) : 
where we can type INNJUTOCH with: 

A hp ATP : A h s . x : (p Tom q> To m Q from 

(TRlocOut) 

A h„ Mit A : A h„ R a \ : {o} • 

£ - ' (TPAR) 

A hp Mit A I Raf ■ {o} 
A h s 7/m : <p rim {o} C (p Tim 

(TCOMP) 

A hp (77m) J (Mit A \Rat)~® 



4 Related Works and Conclusions 



K-calculus is a formalism proposed by Danos and Laneve [11] that idealises protein-protein interactions 
using graphs and graph-rewriting operations. A protein is a node with a fixed number of sites, that may 
be bound or free. Proteins may be assembled into complexes by connecting two-by-two bound sites 
of proteins, thus building connected graphs. Collections of proteins and complexes evolve by means 
of reactions, which may create or remove proteins and bounds between proteins: K"-calculus essentially 
deals with complexations and decomplexations, where complexation is a combination of substances into 
a new substance called complex, and the decomplexation is the reaction inverse to complexation, when a 
complex is dissociated into smaller parts. These rules contain variables and are pattern-based, therefore 
may be applied in different contexts. Even if fc-calculus does not deal with membranes, we have in 
common the use of variables and contexts for rule application. Moreover, both approaches emphasise 
the key rule of the surface components, in proteins (fc-calculus) or membranes (CLSLR), for biological 
modelling. 

P-Systems JT7 1 are a biologically inspired computational model. A P-System is formed by a mem- 
brane structure: each membrane may contain molecules, represented by symbols of an alphabet, other 
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membranes and rules. The rules contained into a membrane can be applied only to the symbols contained 
in the same membrane: these symbols can be modified or moved across membranes. The key feature 
of P-Systems is the maximal parallelism, i.e., in a single evolution step all symbols in all membranes 
evolve in parallel, and every applicable rule is applied as many times as possible. Locality and intrinsic 
parallelism of rules are also present in our approach, but in CLSLR the level of parallelism is not neces- 
sarily maximal, and moreover not only molecules but also rules can be created, deleted or moved across 
membranes. In both approaches the local rules cannot describe some possible biological behaviours such 
as fusion, deletion or creation of membranes. P-Systems are not so flexible in the description of new ac- 
tivities observed on membranes without extending the formalism to model such activities. In CLSLR 
this limit is overcome by global rules, that contain generic patterns. 

In rewrite system models, the term (describing the systems under consideration) and the list of rules 
(describing the system's evolution) could be considered as separate (written on two different sheets of 
paper). In this work, we have presented a calculus with global (separate from the system) and local 
(dynamic and system intrinsic) rewrite rules. While global rules can, as usual, be applied anywhere in a 
given term, local rules can only be applied in the compartment on which they are defined. Local rules 
are equipped with dynamic features: they can be created, moved and erased. 

As it happens for P-Systems, local rules are intrinsically parallel. Indeed, expressing rules that are 
local to well delimited compartments, and with the possibility to define systems with multiple, parallel, 
compartments, naturally leads to the definition of a parallel semantics. 

As a future work, in the lines of (UEKH, we plan to investigate how to adapt this model with a 
quantitative semantics, also studying the limits and constraints imposed by a parallel semantics. 
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Lemma A.l (Inversion Lemma). 1. If Ah s s : (p, then (p = 0. 

2. If A h s x '■ <P> then ^ : 9 S A. 

3. If A h s a : (p, then a : (p € A. 

4. If A h s SP ■ SP' : (p, then A h s SP : (pi, A h s SP' : <?>> and <p = q>\ U (pi. 

5. If A h p X : X, then X:TGA. 

6. IfA\-pSP:X, thenx = %. 

7. If A \- p L\ h-> L2 : T, then % = features(Li H> L2) U %' and A \- p L2 : t '. 

8. If A hp L\ s > ^ Lf 2 : x, then T = {0} • %', Ah p L 2 : %', A h s S\ : <pi, A h, S 2 : <p 2 and (p { C ^ 

9. If A hp L\ Sl 1-^ L,f 2 : T, then % = {i} U r', A h p L 2 : (p ■ X 1 , A h s Si : (p h A h s S 2 : (p 2 and (p U (pi C (p 2 . 

10. If A hp P \ P' : x, then x = X\ U x 2 , A h p P : Xi and A h p P' : x 2 . 

11. If A hp (SP) J P : x, then A h s SP : <p, A h p P : <p' • % and <p' C <p. 

12. If A h gr P\h^P 2 : ok, then Ah p Pi : Ti, A h p P 2 : x 2 and X 2 C Ti. 

Proof. Immediate from the typing rules in Figures |4] and [51 □ 

Lemma A.2. If A h p E[P] : x then 
1. A hp P : Xofor some To, and 
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2. ifP' is such that A \- p P' : x' with x' C T , then A \- p E[P'\ : x" with x" C x. 

Proof. By induction on the definition of contexts. 

• If E = □, then E[P] = P, and so A \- p P : x. Since in this case E[P'\ = P', and A \- p P' : x' with 
x' C T by hypothesis, then A h p £[P'] : t'. 

• If £ = E' | T, then E[P] =E'[P] \ T. From LemmaEO][lO]> we derive A h p E'[P] : Ti and A h p T : T 2 , 
with Ti U t 2 = By induction hypothesis on E'[P] we get A \- p P : To and A h p E'[P'] : tJ with 
t{ C Ti. Applying rule (Tpar) we conclude A \- p E[P'] : x" with x" = x{ U T2, and then t" C t. 

• If £ = (S)°j£", then E[P] = (5)°j£'[P]. From Lemma lAUHI]) we derive Ah p S: <po, and 
A hp E'[P] : (p ■ X, with (p C <p . By induction hypothesis on E'[P] we get A h p P : To, and A \- p 
E'[P'} : <p' • x' with (p' ■ x' Q (p ■ x. Applying rule (Tcomp) we conclude A h p E[P'] : t', with t' Q x. 

□ 

Lemma A.3. If<T£ Za» h s 5Pa : <p //'cotJ o«/j j/A h s : 9. 

Proo/ (<^) By induction on A \- s SP : <p. Consider the last applied rule. 

• If the rule is (tsvar), the proof follows from a € S A . For rules (tseps) and (tselm), the fact that SP 
is a term implies that SPo = SP, and, moreover, it is typable from the empty environment. 

• Rule (tsseq). In this case SP = SPi ■ SP2, and from Lemma lATTBl ) we derive A h $ SP[ : <pi, A h s . 
SP2 '■ <Pi, and (p = (p\ U (p2- By induction hypotheses on SPi and SP2 we get h s SPiO : (pi and 
h, SP2O : <j[>2- Therefore, since 5Pia ■ SP2O = (5Pi • SP2)o, applying the rule (tsseq) we conclude 
h s (SP l -SP 2 )a:( P . 

(=^) By induction on SP. 

• If SP = X, the proof follows from a G £a- If SP = e or SP = a we use weakening. 

• Let SP be SP • SP 2 . Since (SPi • SP 2 )o = SPi a • SP 2 a, from LemmalAZTJU) we derive 9 = <pi U qy>, 
\- s SP\G : <pi, and h s 5P20 : 92. By induction hypotheses we get A h v 5Pi : <pi, and A h s SP2 : 92- 
Applying rule (TSseq) we conclude A h s 5Pi • SP2 : 9. 

□ 

Lemma A.4. /fa G Za> \- p Pa : t jfanrf only if A \- p P : x. 

Proof. (<;=) By induction on A h p P : x. Consider the last applied rule. 

• If the rule is (Tvar), the proof follows from a G Za. For rules (trloc), (trlocOut), (trlocIn) the fact 
that P is a rule implies that Pa = P and, moreover, it is typable from the empty environment. For 
rule (Tseq) if P is a sequence pattern then also Pa is a sequence pattern, and then we can apply rule 
(Tseq) with the empty environment. 

• If the rule is (Tpar), then P = Pi | P 2 , and from Lemma lATTIfTOl ) we derive A h p Pi : Ti , A \- p P 2 : T 2 , 
and T = Ti U t 2 . By induction hypotheses on Pi and P 2 we get \- p P\0 : Ti, and h p P 2 a : T 2 . 
Therefore, since Pia | P 2 a = (Pi | P 2 )a, applying the rule (Tpar) we conclude h p (Pi | P 2 )a : x. 

• If the rule is (Tcomp), then the proof is similar using Lemmas [A. llfTTT ) and |A.3| for the first premise. 
(=^) By induction on P. 

• If P = X, the proof follows from a G £a- If P is a sequence pattern, then also Pa is a sequence 
pattern, and we can apply the rule (Tseq). If P is a rule, then P = Pa. 
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• Let P be P = Pi | P 2 . Since (Pi | P 2 )o = P x aP 2 a, and the fact that h p (Pi | P 2 )a : T, from Lemma 
lA.llfTOl ) we derive \- p P\0 : T\ , h p P 2 o : T 2 , and z = Z\ U T 2 . By induction hypotheses on Pi and P 2 
we get A \- p Pi : Ti and A h p P 2 : T 2 . Applying rule (Tpar) we conclude A h p (Pi | P 2 ) : z. 

• If P = (SP) J P' the proof is similar using Lemmas [A. 1 tTTTT t and [A3] for the first premise. 

□ 

Proof of Theorem [33] (Subject Reduction) 

By cases on the reduction rules. 

Rule (TGR) 

From Definition ED T = E[P y o], T' = E[P 2 o], and a G I A . By hypothesis h p T : z and A h gr 
Pi i— > P 2 : ok. Therefore, Lemma |AT2l fT1) implies h p Pi a : Ti for some Zi, and from Lemma [A~4l we 
derive A \- p Pi : Ti. From A h gr Pi i-> P 2 : ok, Lemma |A~T|[T2] > implies A h p P 2 : T 2 with T 2 C Ti. We 
can apply Lemma IAT41 obtaining h p P 2 a : T 2 . From Lemma lAT2ll2l) we conclude A h p P[P 2 a] : z' 
for some t' C z. 

Rule (TLR) 

From Definition 13. II T = E[L\ (->• L 2 | Ljcr], P' = P[Pi h-> L 2 | P 2 a], and a € E A . By hypothesis 
\- p T : T. Therefore, Lemma lA!2l fT1) implies h p Li i-> P 2 j Pi a : To for some To. Since a € E A , 
Lemma |A!4l implies A h p Li i->- L 2 | L\ : To. From Lemma [ATTUTOt we derive A h p L\ i-> L 2 : Tq, 
and A h p Li : Ti with Tq U Ti = To. By Lemma IATTTP71) we derive A h p P 2 : T 2 with T 2 C Ti. Lemma 
IA.4l implies h p P 2 a : T 2 , then from (Tpar) we derive h p L\ i-> L 2 | L 2 C7 : T3 with T3 = Tq U T 2 . Since 
T3 E To we can apply Lemma IAT2M 2I) obtaining h p P[Pi 1-4 P 2 | P 2 a] : t' with z' C T. 

Rule (lr-out) 

From Definitional] P = P[(Si(t) J (Pf 1 h-> pf 2 | L x o | P )], P' = E[L 2 o | {S 2 o)° \ (Pf 1 1-4 
l\ S2 I P )], and a G E A . By hypothesis h p P : T. Lemma [A2BB implies h p (Sict) C J (Pf 1 1— > 
P| S2 I Pi a I T Q ) : T , and, since a G I A , Lemma[A3]implies A h p (5i)° J (l\ S[ ^ Lf 2 \ L y | P ) : T . 
By Lemma [ATTHTTb we get A h s S\ : <pi, and A h p Pj Sl h4 Pj | Pi | Po : (po • To where <jf>o ^ 9i- By 
Lemma lATTHTOb we have A h p Po : (p ■ Ti , and A h p p{ Sl h-> P^ 2 :<?>'• T 2 for some <p • Ti U q>' ■ z 2 E 
<po • To- Lemma |A. Ilf8l) implies (p' = {o}, A h p P 2 : T 2 and A h s S 2 : (ffy where <pi C <p 2 . Since 
a G £ A , Lemma lA!4l implies h p P 2 a : T 2 , \- s S 2 o : (p 2 , A h p pj 5 ' (->• P^ 52 : {0} • T 2 , and h p Po : <p • Ti . 
Using these premises, we apply rule (Tpar) deriving h p Pj 5 ' i-> P^ 2 | Po : {0} • T 2 U (p ■ Ti, and 
then h p (S 2 a)°J (pf 1 ^ P^ 2 | P ) : Tl u t 2 by rule (Tcomp). Finally h p P 2 a | (S 2 o)° J (pf 1 ^ 
P 2 ^ 2 I Po) : Ti U T 2 by rule (Tpar): since Ti U T 2 C To, we can apply the Lemma lAT2ll2l ). obtaining 
h p P[P 2 a I (S 2 o) } (Lf 1 H4 pf 2 I P )] : z' with z' C t. 

Rule (lr-in) 

From Definitional] P = P[Pf 1 ^Pf 2 | Pi a | (Si<r)° J T ], T' =P[Pf 1 M.pf 2 | (S 2 o)° J (P | P 2 a)], 
and a G £ A - By hypothesis h p P : t. Lemma[A2]Q) implies h p p{ 5 ' h-> P^ 2 | Pi a \ (Si a) J P : T , 
and, since a G I A , Lemma [A3] implies A h p pj Sl pf 2 | Pj j (Si)°J P : T . By Lemma |A~T1[T0] ) 
we have A h p Pj Sl 1-4 P^ 2 : Ti and A h p (Si) J Po : T 2 for some Ti U t 2 C t<). Lemma |A~T] |9]> implies 
Ti = {i} U T 3 , A h s 5i : <pi , A h. s 5 2 : (jf^, and A h p P 2 : <p • T 3 , where <p U q>\ C <p 2 . By Lemma lATTIfTTT ) 
A h p Po : (p' ■ z 2 for some (p' Q (pi. Since a G E A , Lemma lA~4l implies h s 5iC7 : (pi, h v 5 2 a : (p 2 , 
\- p L 2 o : <p • T3, h p p{ Sl (->■ P^ 2 : {i} U T3, and h p P : <p' • T 2 . Using these premises, we apply rule 
(Tpar) deriving h p P 2 a | Po : <p • T3 U <p' • T 2 , and then h p (S 2 o)° \L 2 a \ Po : T3 U T 2 by rule (Tcomp). 
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Finally \- p L\ 1 h-> Lf 1 \ (S 2 o) \ L 2 G \ T : Ti U z 2 , because X\ = {i} U T 3 , by rule (Tpar): since 
Ti U T 2 E T , we can apply the Lemma lAT2lf2T > obtaining h p E[l\ Si i-> L^ 2 [ (^a) J L 2 a | 7b] : t' 
for some t' C t. 



